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Abstract. We present thermal expansion and longitudinal magnetostriction data 
for cubic RA1 3 (R = Tm, Yb, Lu) single crystals. The thermal expansion coefficient 
for YbAl 3 is consistent with an intermediate valence of the Yb ion, whereas the 
data for TmAl3 show crystal electric field contributions and have strong magnetic 
field dependencies. de Haas-van Alphen-like oscillations were observed in the 
magnetostriction data of YbAl 3 and LuAl 3 , several new extreme orbits were measured 
and their effective masses were estimated. Zero and 140 kOe specific heat data taken on 
both LuAl 3 and TmAl 3 for T < 200 K allow for the determination of a CEF splitting 
scheme for TmAl 3 . 
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1. Introduction 

Rare earth trialuminides, RAI3, have been studied for several decades. The crystal 
structure of these materials is very sensitive to the rare earth ionic radius and it changes 
from hexagonal for light rare earths to cubic (CU3AU - type) for TH1AI3, YbAl3, LUAI3 
(and SCAI3) PQ, with ErAl3 reported to have been synthesized in both crystallographic 
variants [1, 2, 3J. TmAl 3 was reported to have a singlet ground state [31 H]: the magnetic 
susceptibility follows the Curie- Weiss law with the magnetic moment corresponding to 
that of Tm 3+ at high temperatures and becomes a temperature-independent Van Vleck 
type at low temperatures. YbAl3, an intermediate valence compound [5j with the high 
Kondo temperature, T K , of 600-700 K, recently experienced a revival of interest to 
its physical properties [HI [7J [H [9] based, in part, on the uncovering of a second, low 
temperature, energy scale, Fermi liquid coherence (T co h « 40 K), and slow crossover 
between T K and T coh [TO] . 

Whereas many physical properties were measured and analyzed for cubic trialu- 
minides, thermal expansion and magnetostriction data for these materials, in particular 
at low temperatures, appear to be absent from the literature (except for an early report 
[TTj on the thermal expansion of LUAI3 and YbAl3 between ~ 90 K and ~ 800 K). In 
addition, for TmAl 3 the exact crystal electric field (CEF) splitting of the Hund's rule, 
ground state multiplet has remained unresolved [U [12], [131 E]> i n P ar t due to a lack 
of high temperature and high magnetic field specific heat data. In this work we report 
extensive measurements of specific heat as well as thermal expansion and magnetostric- 
tion data for RA1 3 (R = Tm, Yb, Lu) single crystals so as to compare with other salient 
physical properties and to inquire into the effects of the CEF and intermediate valence 
(IV) on these characteristics. Simple crystal structure and the ability to synthesize 
large, high quality, singly crystals facilitate this endeavor. 

2. Experimental methods 

Large (> 0.1 cm 3 ) single crystals (Fig. [lj inset) of RAI3 (R = Tm, Yb, Lu) and 
R'o.iLuo.9Al 3 (R' = Er, Tm, Yb) were grown from aluminum-rich binary (R - Al) or 
ternary (R' - Lu - Al) melts using a self-flux method [15]. High purity R (Ames Lab) 
and Al were placed, in atomic ratios of R0.12AI0.88 or (Luo.gRo.Oo.i^Alo.ss; in alumina 
crucibles. These were sealed in fused quartz ampoules under a 1/3 atmosphere partial 
pressure of Ar. The ampoules were heated to 1050°C and then cooled to 675°C over 
150 - 160 hours, at which point they were removed from the furnace and the excess Al 
was decanted. For R'o.iLuo.gALj with R' = Er, Tm the nominal concentration of R' was 
corroborated by the Curie- Weiss fit of the high temperature susceptibility. The inset to 
Fig. [1] shows crystals with clear, cubic morphology and mirrored (100) facets. Linear 
dimensions as large as 7 mm were commonly achieved. 

Thermal expansion and magnetostriction were measured using a capacitive 
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dilatometer constructed of OFHC copper; a detailed description of the dilatometer is 
presented elsewhere [16]. The dilatometer was mounted in a Quantum Design PPMS-14 
instrument and was operated over a temperature range of 1.8 to 300 K in an applied 
magnetic field up to 140 kOe. The samples were mounted in such a way that thermal 
expansion was measured along the [100] direction. The applied magnetic field was also 
along [100] such that H \\L\\ [100]. The crystals were cut and polished so that the typical 
distance between the parallel (100) surfaces of the samples was L[ioo] ~ 2 — 3 mm. Heat 
capacity of the samples was measured using a hybrid adiabatic relaxation technique of 
the heat capacity option in a Quantum Design PPMS-14. 

Thermodynamic properties of materials are frequently analyzed using the concept 
of a Griineisen function (or a Griineisen parameter) [17]. For a single salient energy 
scale, e, the Griineisen parameter, 7, is defined as 7 = — dine /din V, where V is a 
molar volume. Using thermodynamic relations, we can obtain 7(T, V) = (3V/xsC p , 
where (3 is a volume thermal expansion coefficient (f3 = (din V/dT)p), xs is an 
adiabatic compressibility (xs — —(din V/dP)s) and C p is a heat capacity at a constant 
pressure. For cubic materials 7(T, V) = 3aV/x,sC p , where a is a linear thermal 
expansion coefficient (a = (din a/dT)p, a - is a lattice parameter). Sometimes, in 
the analysis of experimental data, lacking the temperature-dependent compressibility 
data, the temperature dependence of the Griineisen parameter can be approximated 
[T9] as being proportional to (3/C p (or 3a /C p for cubic materials) under the assumption 
that the relative temperature dependence of xs is significantly smaller then that of 
thermal expansion coefficient or heat capacity. We will adopt such approach in this 
work. 

If more than one contribution to the thermodynamic properties is present (e.g. vi- 
brational, electronic, magnetic, etc.), the Griineisen parameters are not additive, rather 
the Griineisen parameter for the material is an average, weighted by the heat capacity 
contribution of each component [IT]: 7 = ^7 r C r / ^ C r . 

r r 

3. Results and discussion 

3.1. Thermal expansion and heat capacity 

Zero field, temperature-dependent, linear (L||[100]) thermal expansion coefficients for 
RAI3 (R = Tm, Yb, Lu) (together with the literature data [H] of polycrystalline Cu 
for comparison), are shown in Fig. [TJ In the overlapping temperature regions (T > 90 
K), the thermal expansion coefficient values for LUAI3 and YbAl3 deduced from the 
the lattice parameters as a function of temperature data [11] are comparable to our 
results. At room temperature, thermal expansion values of the non-hybridizing TmAl3 
and LUAI3 are very similar. The differences in a(T) between these two materials, on 
cooling, and in particular a peak in a(T) of TmAl 3 at ~ 20 K, are probably related 
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to the CEF contributions to the thermal expansion of TmAl 3 . a(T) of YbAl 3 is lower 
than that of its non-magnetic analogue, LuAl 3 , over the whole temperature range. Such 
behavior is consistent [20J with YbAl 3 being an Yb-based, intermediate valence mate- 
rial with a high, Tk 3> 300 K, Kondo temperature. Qualitatively, such behavior can be 
understood by noting that the fractional Yb valence in YbAl 3 increases with increase 
in temperature (in the temperature range of our measurements) [21] and that the ionic 
radius of Yb 3+ is smaller than that of Yb 2+ . 

Figure [2] presents the temperature-dependent linear thermal expansion coefficient 
(a), the heat capacity (C p ) and the ratio a/C p for LuAl 3 . a [100] CO and C P (T) have 
similar temperature dependencies. The ratio of these two quantities (that, as men- 
tioned above, is likely to approximate the temperature dependence of the Griineisen 
parameter) is practically constant down to ~ 50 K, only rising at lower temperatures 
and manifesting a peak at ~ 8 K. It is not unusual to observe a temperature-dependent 
Griineisen parameter, even in simple, non-magnetic, metals [XT] . This is due to the dif- 
ferent temperature dependencies of the electron and phonon contributions to the ther- 
modynamic properties, with these contributions becoming comparable in magnitude at 
low temperatures. Additionally, the error bars in a/C p may be somewhat enhanced at 
low temperatures as a result of calculating the ratio of two, diminishingly small numbers. 

For TmAl 3 (Fig. [3]) both thermal expansion and heat capacity exhibit a broad peak 
in the 15 — 20 K temperature range, apparently related to CEF effects. This peak is still 
observable in the a/C p curve, however, if the data are plotted as Aa/AC p , (i.e. after 
subtraction of the non-magnetic background) it is not seen, confirming the same origin 
of the feature in C p and a. The Aa/AC p shows a close to linear, slowly changing with 
temperature, behavior, probably representing rather well the behavior of the magnetic 
contribution to the Griineisen parameter. 

The peak in heat capacity evolves with applied field (Fig. [4]) reflecting the lifting 
of the degeneracy of TmAl 3 energy levels. A variety of differing CEF schemes, based 
on fits of different experimental quantities (C P (T), or x(T), or high field M(H)) were 
suggested for TmAl 3 in the literature [4[ [T2J [13j [14]. To analyze our magnetization, 
susceptibility and specific heat data we will use the approach delineated in Ref. [22J. 
Since Tm 3+ ions in TmAl 3 have octahedron type of coordination, both parameters, W 
and x, of the CEF hamiltonian for cubic symmetry (see Ref. [22] for definitions of these 
parameters and detailed discussion) are negative and the ground state is T 2 (set 1) or 
Ti (set 2). Temperature dependent susceptibility and magnetization isotherms up to 70 
kOe (Fig. [5]) allow for similar quality fits for CEF schemes with either ground states. 
Heat capacity data are better fit with the CEF scheme corresponding to T 2 ground 
state (Fig. [6]) (for both sets the W and x values from fits of temperature-dependent 
susceptibility and magnetic isotherms were used). The W and x values for both sets 
are listed on the graph. An applied magnetic field changes the CEF splitting via the 
Zeeman term. The H = CEF levels plus this Zeeman term describe H = 140 kOe 
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data well (Fig. EJ. 

Fig. [7] shows the simulation of the heat capacity behavior using the various CEF 
splitting schemes proposed in the literature jU [121 EE EH] including two of the three 
schemes [I] presented as indistinguishable given their low temperature data. Set 2 in 
Ref. [4] and the parameters from Ref. p2] are close to our results (set 1) and describe 
the heat capacity data reasonably well, whereas the parameters from Refs. [T3"l [14"] and 
set 1 in Ref. [4] fail as reasonable approximations of the experimental data. 

An applied field changes the behavior of the thermal expansion coefficient of TmAl 3 
dramatically (Fig. [8]). In between 25 kOe and 50 kOe the low temperature maximum 
turns into a minimum. This minimum deepens in higher fields and reaches ~ —1.7- 10~ 5 
K _1 near 15 K for H = 140 kOe. Clearly, as a result of the changing CEF energy scales, 
the Griineisen parameter will be significantly different in the applied field. Qualitatively 
similar, but approximately an order of magnitude lower in size, field-induced changes 
in aq 10 o](T) are observed in Tmo.1Luo.9Al3 (cf. two insets to Fig. EJ). These perceptible 
field dependences of a[i 00 ](T) in both TmAl 3 and Tm 01 Luo.9Al 3 are consistent with a 
CEF- related, single-ion effect. 

To illustrate the complexity of the effect of magnetic field on the thermal expan- 
sion coefficient, data for Er .iLuo.gAl 3 is shown in Fig. [91 As opposed to Tmo.1Luo.9Al3 
(Fig. El inset) a positive, broad feature in a[ W0 ](T) below ~ 40 K grows with increasing 
applied magnetic field. Transverse thermal expansion measurements in an applied field 
and some knowledge the CEF scheme for Ero.1Luo.9Al3 would be needed for comparative 
analysis of the effects of magnetic field on these two different materials. 

For YbAl 3 , a [100] CO and C P (T) (Fig. [10]) have similar temperature dependencies, 
without any apparent striking features. The ratio of these two quantities is practically 
constant down to ~ 30 K and then decreases at lower temperatures. Aa/AC p is linear 
(and close to constant) at higher temperatures and then, on further cooling, decreases 
and passes through broad minimum. Given that the Kondo temperature is above room 
temperature (5], it is tempting to try to connect these changes in Aa/AC p to the 
emerging second, low temperature (~ 30 — 40 K), scale in YbAl 3 [10], however at 
this point there is no clear evidence for such conjecture and more studies are required. 
Additionally, we cannot exclude this low temperature behavior to be the result of an 
interplay between phonons and (enhanced) electronic degrees of freedom. Specific heat 
data show an elevated electronic specific heat coefficient, 7 ~ 50 mJ/mol K 2 , for YbAl 3 , 
significantly higher than that for LuAl 3 , in consistence with the intermediate valence 
nature of YbAl 3 and similar to recently reported data [23] . 

a [ioo](T) for YbAl 3 is not significantly affected by an applied magnetic field of 140 
kOe (Fig. [TT1) . This is not surprising for a material with a high, 600 — 700 K, Kondo 
temperature. At low temperatures (T < 20 K) YbAl 3 exhibits a region of negative 
thermal expansion, «[ioo](T) < 0. For Ybo.iLuo.gAl 3 the region of negative thermal 
expansion does not appear (above 1.8 K) since a[i 00 ](T) is apparently dominated by the 
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contribution from the LuAl 3 matrix. Both, in pure YbAl 3 and diluted Ybo.1Luo.9Al3, 
the temperature dependent thermal expansion is below the one for the non-magnetic 
analogue, LUAI3, consistent with an intermediate valence character of the Yb ion in pure 
and diluted material. 

3.2. Magnetostriction and quantum oscillations 

In LuAhj the base temperature (1.8 K) magnetostriction is rather small, however, start- 
ing at fields below 40 kOe de Haas-van Alphen (dHvA) like oscillations in magnetostric- 
tion (Fig. [T2|) are clearly observed. A fast Fourier transform performed on these data 
in the form of AL/L vs. 1/H. Seven dHvA frequencies ranging from ~ 1 MG to 
~ 50 MG were observed (Fig. [T2], inset). The occurrence of quantum oscillations in 
magnetostriction is a known phenomenon [24J, however observations of such oscillations 
are rather rare, since both large, high quality single crystals and sensitive dilatometers 
are required. dHvA oscillations were observed in LUAI3 via magnetostriction up to tem- 
peratures as high as 20 K (Fig. [13]). The temperature dependence of the amplitude of 
these oscillations can be used to evaluate the effective masses of the quasiparticles on the 
corresponding extremal orbits using standard Lifshitz - Kosevich formula [25, 26, 27J. 
The obtained dHvA frequencies and corresponding effective masses, m*/m , together 
with the literature data [28] obtained by a conventional magnetic susceptibility field- 
modulation technique are shown in Fig. [TH Of the seven orbits observed in this work, 
four (above 10 MG) are consistent with experimental or theoretical literature data and 
three frequencies (F < 10 MG, m* /m$ < 0.5) are new. 

The magnetostriction of YbAl 3 is also rather small: AL/L 3 • 10~ 9 at T = 1.8 
K, H = 140 kOe. In high fields dHvA-like oscillations are also observed (Fig. [T5l) . The 
amplitude of the oscillations is significantly smaller (few orders of magnitude for the 
leading frequency) than that for LuAl 3 and they reduce to the level of noise above ~ 15 
K. Two frequencies and their effective masses were identified from our measurements 
and the results are plotted together with the literature data in Fig. [161 It is worth 
noting that, for the orbits detected by magnetostriction the electronic masses in YbAl 3 
are similar to those found for LUAI3. Significant mass enhancement appears to occur 
only for higher frequencies. 

Longitudinal magnetostriction of TmAl 3 is orders of magnitude higher than that 
of YbAl 3 and LuAl 3 : AL/L « 7 ■ 10" 4 at T = 1.8 K, H = 140 kOe (Fig. HTJ). In 
the intermediate and high fields the magnetostriction is approximately proportional to 
H 2 (Fig. [T71 inset). Detailed analysis of the magnetostriction [3TJ [32j [33] requires 
knowledge of the details of the Tm3 + ion-lattice interactions. 

We can not resolve quantum oscillations of magnetostriction in TmAl 3 , although 
there are no indications of inferior crystal quality. We need to mention though that the 
much higher, monotonic magnetostriction background makes observation of quantum 
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oscillations extremely difficult. Using a traditional magnetic susceptibility modulation 
technique seven fundamental dHvA frequencies were observed in TmAl 3 for £f||[100] 
|34j . these frequencies being similar to the ones found for LUAI3 and YbAl3. 

4. Summary 

Temperature-dependent thermal expansion coefficients were measured between 1.87 K 
and 300 K for RA1 3 (R = Lu, Yb, Tm) along [100] direction. In YbAl 3 the intermediate 
valence of the Yb-ions results in a (T) for this material being consistently lower than a (T) 
of the non-magnetic analogue, LuAl 3 , with a negative thermal expansion region at low 
temperatures. Thermal expansion coefficient, a(T), of TmAl 3 manifests contribution 
from the CEF effects and, as opposed to YbAl 3 and LuAl 3 , displays strong field 
dependence below ~ 100 K. 

Magnetostriction of YbAl 3 and LuAl 3 [100]) shows clear dHvA-like 

oscillations up to temperatures as high as above 20 K for LuAl 3 . Several new dHvA 
frequencies were measured (three for LuAl 3 and two for YbAl 3 ) and their effective masses 
were estimated. For these orbits the electronic masses in YbAl 3 are similar to those found 
for LuAl 3 . TmAl 3 shows very large, exceeding that of YbAl 3 and LuAl 3 by several orders 
of magnitude, magnetostriction at low temperatures. 

In addition, a CEF splitting scheme for TmAl 3 with the T 2 level as the ground 
state is strongly suggested based on the analysis of the heat capacity to significantly 
higher fields and temperatures. This result removes a two decade ambiguity as to the 
CEF level scheme for this material. 
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Figure 1. (Color online) Temperature-dependent linear thermal expansion coefficients 
along [100] direction of RAI3 (R = Tm, Yb, Lu) in zero applied field. Data for 
polycrystalline Cu [TH] are shown as dashed line for comparison. Inset: LUAI3 crystal 
over a mm scale. 
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Figure 2. (Color online) Temperature-dependent linear thermal expansion coefficient, 
heat capacity and the ratio a/C p for LuAl 3 in zero applied field. 
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Figure 3. (Color online) Temperature-dependent linear thermal expansion coefficient, 
heat capacity and the ratio a/C p for TmAl3 in zero applied field. Additionally, on the 
lower panel Aa/AC p = (a(TmAl 3 ) - a(LuAl 3 ))/(C p (TmAl 3 ) - C p (LuAl 3 )) is plotted. 
Dashed line is a guide for the eye. 
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Figure 4. (Color online) The low temperature part of heat capacity for TmAi3 in 
applied fields H = 0, 50, 100, 140 kOe. Field was applied along [100] axis. Inset: data 
for H = and 140 kOe up to 200 K. 
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Figure 5. (Color online) (a) Temperature-dependent magnetic susceptibility, M/H, 
inverse magnetic susceptibility, H/M, taken at H — 1 kOe, -ff||[100] and (b) 
magnetization isotherms, M(H), for T — 1.85 K and 300 K for TmAl3. Curie 
- Weiss fit of inverse magnetic susceptibility and CEF simulations of H/M and 
M(H) for two possible sets of parameters (set 1: W = —1.011, a; = —0.298 set 2: 
W = —2.125, x = —0.682) are show as lines. 
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Figure 6. (Color online) Magnetic contribution to the heat capacity of TmAi3, 
C m = C^ mAh - C^ uAh , in zero and 140 kOe (F||[100]) applied field. CEF simulations 
for two possible sets of parameters are show as lines. Corresponding W and x values 
are listed for both sets are listed. CEF scheme for set 1 is illustrated in the right side 
of the plot. 
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Figure 7. (Color online) Magnetic contribution to the heat capacity of TmAi3, 
C m = C^ mAh -C^ uM \ in zero and 140 kOe (F||[100]) applied field. CEF simulations 
using set 1 from this work and the literature data are show as lines. The following 
W and x values (listed as {W, x}) were used: this work, set 1: { — 1.011, —0.298}; Rcf. 
g], set 1: {1.95,0.82}, set 2: {-0.993,-0.282}; Ref. [12j: {-0.94,-0.27}; Ref [B]: 
{3.24, -0.8}; Ref [H]: {-2.316, -0.827} 
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Figure 8. (Color online) Temperature-dependent linear thermal expansion coefficient 
for TmAl 3 in applied fields of 0, 25, 50, 75, 100, 125, 140 kOe. Insets: low temperature 
part of a(T) in 0, 50, 100, 140 kOe applied fields for TmAl 3 and Tmo.1Luo.9Al3. 
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Figure 9. (Color online) Temperature-dependent linear thermal expansion coefficient 
for Ero.1Luo.9Al3 in applied fields of 0, 50, 100, 120, 140 kOe. Insets: low temperature 
part of a(T). 
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Figure 10. (Color online) Temperature-dependent linear thermal expansion 
coefficient, heat capacity and the ratio a/C p for YbAl 3 in zero applied field. 
Additionally, on the lower panel Aa/AC p = (a(YbAl 3 ) - a(LuAl 3 ))/(C p (YbAl 3 ) - 
C p (LuAl 3 )) is plotted. Inset: low temperature C p /T vs. T 2 for LuAl 3 and YbAl 3 . 
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Figure 11. (Color online) Low temperature part of linear thermal expansion 
coefficient for YbAl 3 and Ybo.1Luo.9Al3 in applied fields of and 140 kOe. For 
comparison, data for LUAI3 are shown as dashed line. 
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Figure 12. High field part of longitudinal (H \\L\\ [100]) magnetostriction in LuAi3 
measured at 1.8 K. Inset: fast Fourier transform of the corresponding AL/Lq vs. l/H 
data. Arrows mark the de Haas-van Alphen frequencies. 
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Figure 13. (Color online) Longitudinal magnetostriction in LUAI3 for if||[100] at 
different temperatures, from 1.8 K to 20 K plotted as a function of 1/H. Constant 
background was subtracted from the data. 
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Figure 14. (Color online) dHvA frequencies and corresponding effective masses for 
LuAl 3 , i?||[100]. Filled circles - this work, other symbols - literature data [55]: open 
circles - experiment, triangles - theory. Note that Sakamoto et al. experimentally 
identified several other dHvA frequencies but did not determine the corresponding 
effective masses. These experimental frequencies are shown by arrows. 
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Figure 15. Longitudinal magnetostriction in YbAl 3 for i?||[100] at T — 1.8 K plotted 
as a function of 1/H. Constant background was subtracted from the data. 
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Figure 16. (Color online) dHvA frequencies and corresponding effective masses for 
YbAl 3 , H\\ [100]. Filled circles - this work, other symbols - literature data: open circles 
- experiment [29], triangles - theory [30]. Data for LUAI3 (this work) are shown as 
asterisks for comparison. 
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Figure 17. (Color online) Longitudinal magnetostriction in TmAi3 for ff||[100] at 
different temperatures, from 1.8 K to 100 K. Note that within the scale of the plot, 
1.8 K and 5 K data coincide. Inset: the same data on a log - log plot. Dashed lines 
correspond to AL/L cx H 2 functional behavior. 



